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COSMIC EVOLUTION OF STELLAR DISK TRUNCATIONS: 
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ABSTRACT 

We have conducted the largest systematic search so far for stellar disk truncations in disk-like 
galaxies at intermediate redshift (z<l.l), using the Great Observatories Origins Deep Survey South 
(GOODS-S) data from the Hubble Space Telescope - ACS. Focusing on Type II galaxies (i.e. down- 
bending profiles) we explore whether the position of the break in the rest-frame B-band radial surface 
brightness profile (a direct estimator of the extent of the disk where most of the massive star formation 
is taking place), evolves with time. The number of galaxies under analysis (238 of a total of 505) is an 
order of magnitude larger than in previous studies. For the first time, we probe the evolution of the 
break radius for a given stellar mass (a parameter well suited to address evolutionary studies). Our 
results suggest that, for a given stellar mass, the radial position of the break has increased with cos- 
mic time by a factor 1.3±0.1 between z~l and z~0. This is in agreement with a moderate inside-out 
growth of the disk galaxies in the last ~ 8 Gyr. In the same period of time, the surface brightness 
level in the rest-frame -B-band at which the break takes place has increased by 3.3±0.2 mag/arcsec 2 (a 
decrease in brightness by a factor of 20.9±4.2). We have explored the distribution of the scale lengths 
of the disks in the region inside the break, and how this parameter relates to the break radius. We 
also present results of the statistical analysis of profiles of artificial galaxies, to assess the reliability 
of our results. 

Subject headings: galaxies: structure — galaxies: evolution — galaxies : high redshift — galaxies: 
fundamental parameters 



1. INTRODUCTION 

Early studies of the disks of spiral galaxies (Patter- 
son 1940; de Vaucoulcurs 1959; Freeman 1970) showed 
that this component generally follows an exponential ra- 
dial surface-brightness profile, with a certain scale length, 
usually taken as the characteristic size of the disk. Free- 
man (1970) pointed out, though, that not all disks follow 
this simple exponential law. 

Another repeatedly reported feature of disks is that 
of a truncation of the stellar population at large radii, 
typically 2-4 exponential scale lengths (see e.g. the re- 
view by Pohlen et al. 2004). Van der Kruit (1979) 
and van der Kruit & Searle (1981a, b) first drew atten- 
tion to this phenomenon, which they inferred primarily 
from the major axis profiles of edge-on, late-type spirals. 
Though the term "truncation" is used, not even in the 
original studies was a complete absence of stars beyond 
the truncation radius suggested. In fact, Pohlen et al. 
(2002) showed that the truncation actually adopts the 
form of a fairly sharp change in slope, from the shallow 
exponential of the main disk to a steeper exponential at 
larger radii (see also de Grijs, Kregel & Wesson 2001). 
Erwin et al. (2005) denoted galaxies with this feature 
as Type II objects, generalizing a classification scheme 
initially proposed by Freeman (1970). 

Though the truncation phenomenon appears to be very 
widespread (see for example Pohlen & Trujillo 2006, 
hereinafter cited as PT06), there are many cases in which 
this just does not appear to happen even at extremely 
faint surface brightness levels. To give an example, 
Bland-Hawthorn et al. (2005) found (using star counts) 
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a galaxy (NGC 300) for which the single exponential de- 
cline simply continues down to ~ 10 radial scale lengths 
(~ 30.5 mag/arcsec 2 in r' band). Together with earlier 
measurements by e.g. Barton & Thompson (1997) or 
Weiner et al. (2001) (using surface photometry) this 
provides evidence that indeed there are prototypical ex- 
ponential disks (Type I objects, following nomenclature 
in Freeman 1970). More strikingly, perhaps, there even 
exists evidence for a third type of profile presented by 
Erwin et al. (2005, 2007, 2008) for early-type disks, 
by Hunter & Elmegreen (2006) and PT06 for late-type 
disks ; and for extreme late-type spirals by Matthews & 
Gallagher (1997). In this class, named "Type III" (also 
"antitruncations" ) by Erwin et al. (2005), the inner pro- 
file is a relatively steep exponential, which gives way to 
a shallower profile, beyond the break radius. This profile 
is thus something like the "inverse" of a Type II profile, 
bending "up" instead of "down" beyond the break radius 
( "upbending" profile) . 

In this paper we will concentrate on the Type II galax- 
ies. Several possible break-forming mechanisms have 
been investigated to explain their truncations. There 
have been ideas based on maximum angular momentum 
distribution: van der Kruit (1987) proposed that angular 
momentum conservation in a collapsing, uniformly rotat- 
ing cloud naturally gives rise to disk breaks at roughly 
4.5 scale radii. Van den Bosch (2001) suggested that 
the breaks are due to angular momentum cut-offs of 
the cooled gas. On the other hand, breaks have also 
been attributed to a threshold for star formation (SF), 
due to changes in the gas density (Kennicutt 1989), 
or to an absence of equilibrium in the cool Interstellar 
Medium phase (Elmegreen & Parravano 1994; Schaye 
2004). Following this, and using a semi-analytic model, 
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Elmegreen & Hunter (2006) demonstrated that a double- 
exponential profile may arise from a multi-component 
star formation prescription. The above two "simple" 
scenarios (angular momentum vs. star formation thresh- 
old), however, are challenged by observational results. It 
is becoming more and more clear that galaxies have a 
significant density of stars beyond the break radius. In 
addition, the existence of extended UV disks (Thilker et 
al. 2005; Gil de Paz et al. 2005) and the lack of a clear 
correlation between Ha cut-offs and optical disk breaks 
(Pohlcn et al. 2004; Hunter & Elmegreen 2006) further 
complicates the picture. Even though a sharp star for- 
mation or angular momentum cut-off may explain a disk 
truncation, it does not provide a compelling explanation 
for extended outer exponential components. More elab- 
orate models such as that by Debattista et al. (2006) 
demonstrated that the redistribution of angular momen- 
tum by spirals during bar formation also produces real- 
istic breaks, using collisionless N-body simulations. In a 
further elaboration of this idea, Roskar et al. (2008) have 
performed high resolution simulations of the formation of 
a galaxy embedded in a dark matter halo (see also Bour- 
naud, Elmegreen & Elmegreen 2007; Foyle, Courteau 
& Thacker 2008). They are able to reproduce Type II 
profiles, which they claim to be in good agreement with 
observations. In their model, breaks are the result of 
the interplay between a radial star formation cut-off and 
redistribution of stellar mass by secular processes. 

Independently of the formation mechanism of the trun- 
cation, it seems reasonable to assume that the struc- 
tural properties of the faintest regions of galactic disks 
must be intimately linked to the processes involved in 
the growing and shaping of galaxies. These outer edges 
are easily affected by interactions with other galaxies 
and, consequently, their characteristics must be closely 
connected with the evolutionary paths followed by the 
galaxies (Pohlen & Trujillo 2006; Erwin et al. 2007). 
Together with their stellar halos, the study of the outer 
edges allows the exploration of the "fossil" evidence im- 
printed by the galaxy formation process (de Jong et al. 
2007). 

Furthermore, addressing the question of how the ra- 
dial truncation evolves with z is strongly linked to our 
understanding of how the galactic disks grow and where 
star formation takes place. Perez (2004) showed that 
it is possible to detect stellar truncations even out to 
z ~ 1 . Using the radial position of the truncation (here- 
after Rsr) as a direct estimator of the size of the stel- 
lar disk, Trujillo & Pohlen (2005) inferred a moderate 
(~25%) inside-out growth of disk galaxies since z ~ 1. 
An important point, however was missing in the previ- 
ous analyses: the evolution with redshift of the radial 
position of the break at a given stellar mass. The stellar 
mass is a much better parameter to explore the growth 
of galaxies, since the luminosity evolution of the stellar 
populations can mimic a size evolution (Trujillo et al. 
2004, 2006). Addressing this point here, one aim of the 
present work is to probe whether the galaxies are growing 
from the inside outwards, with star formation propagat- 
ing radially outward with time. Another point to men- 
tion is that while Trujillo & Pohlcn (2005), worked with 
imaging data from the Ultra Deep Field (Beckwith et 
al. 2006), we make use of GOODS-HST/ACS data (Gi- 
avalisco et al. 2004). This field provides a much wider 



sky coverage, and so increases the size of the sample un- 
der study by roughly an order of magnitude, though at 
the cost of using images with somewhat lower signal to 
noise ratio. We take care to tackle the implications of 
this caveat in the paper as well. 

In this work (to which we will refer hereinafter as 
ATB08) we also broach other issues relative to the de- 
tailed characterization of the surface brightness profiles 
(ti-i profiles) . We present results on the distribution of 
the -B-band rest-frame surface brightness at the position 
of the break (hereafter fj,Br)- We discuss the distribution 
of the ratio Rbt / hi , where h\ is the scale length of the 
disk in the region inner to the break. Also included are 
results on the relation between h\ and absolute magni- 
tude and stellar mass. Finally, for the first time in this 
kind of studies, we make extended analysis of the profiles 
of simulated galaxies, in an attempt to better limit the 
accuracy of the reported figures. 

The outline of this paper is as follows. In section 2 we 
describe the galaxy sample and give a brief description 
of the imaging data we use. In Section 3 we describe how 
we obtained the radial surface-brightness profiles of the 
galaxies, and the methodology of the analysis we apply 
for their characterization. Section 4 is focused on the 
completeness of the samples under study and assessing 
the accuracy level of the parameters retrieved. Section 5 
presents results on classification and characterization of 
the /i-r profiles, and more extensive analysis of how the 
Break Radii of Type II disks relate to various other prop- 
erties of the galaxies at different redshifts. In section 6 
we summarize our results and conclusions. Finally we 
include an Appendix on the description and analysis of 
results from an exercise of classification of profiles per- 
formed on a sample of artificial model galaxies. Through- 
out, we assume a flat A-dominated cosmology (Om = 
0.30, O\=0.70, and H =70 km s" 1 Mp C - v ). 

2. DATA & SAMPLE SELECTION 

We have worked with objects present in imaging data 
from HST-ACS observations of the GOODS-South field 
(Giavalisco et al. 2004), which subtends an area in the 
sky of roughly 170 square arcmin. These data are pub- 
licly available 2 , and consist of imaging in the F435W, 
F606TU, F775W and F850LP HST pass-bands, here- 
after referred to as -B435 , Vaoe , 2775 and z 850 , and whose 
total exposure times are approximately 7200, 3040, 3040, 
and 6280 seconds, respectively. The images have an an- 
gular scale of 0.03" /pixel and the FWHM of the PSF in 
z 850 measures 0.09". 

To make the selection of our sample we have profited 
from the fact that the GOODS-South field is within the 
Galaxy Evolution from Morphologies and SEDs imag- 
ing survey (GEMS; Rix et al. 2004). In the red- 
shift range 0.1<z<l.l, GEMS provides morphologies and 
structural parameters for nearly 10,000 galaxies (Barden 
et al. 2005; Mcintosh et al. 2005). For many of these 
objects there also exist photometric redshift and lumi- 
nosity estimates, and SEDs from COMBO-17 (Classify- 
ing Objects by Medium-Band Observation in 17 filters; 
Wolf et al. 2001, 2003). The COMBO-17 team made 
this information publicly available through a catalog with 
precise redshift estimates (with errors <5z/(l+z)~0.02) 

2 http://www.stsci.edu/ftp/science/goods/ 
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for approximately 9000 galaxies down to to# < 24 (Wolf 
et al. 2004). In the same data release were included 
rest-frame absolute magnitudes and colors (accurate to 
~ 0.1 mag). We have also used the stellar mass estimates 
published in Barden et al. (2005), which arc taken from 
Borch (2004), and are deduced from the COMBO-17 
photometric data. 

Barden et al. (2005) conducted the morphological 
analysis of the late-type galaxies in the GEMS field by 
fitting Sersic r 1 /™ (Sersic 1968) profiles to the surface 
brightness distributions. Ravindranath et al. (2004) 
showed that using the Sersic index n as the criterion, 
it is feasible to distinguish between late- and early-type 
galaxies at intermediate redshifts. Late-types (Sab-Sdm) 
are defined as having n<2-2.5. Moreover, the morpho- 
logical analysis conducted by Barden et al. provides the 
information about the inclination of the galaxies. This is 
particularly important, since we want to study the trun- 
cation of the stellar disks in objects with low inclination. 
The edge-on view facilitates the discovery of truncations 
but introduces severe quantitative problems caused by 
the effects of dust and line-of-sight integration which we 
want to avoid (Pohlen et al. 2002). 

We selected objects from the Barden et al. sample 
within the subsequent ranges of parameters: Sersic in- 
dex n<2.5 to isolate disk-dominated galaxies (Barden 
et al. 2005; Shen et al. 2003; Ravindranath et al. 
2004); axial ratio q > 0.5 to select objects with inclina- 
tion < 60°; and Mb < -18.5 magnitudes, as in Trujillo 
& Pohlen (2005). Moreover, only objects with z<l.l 
were selected in order to maintain our analysis in the op- 
tical rest-frame bands. Finally the resulting sample was 
matched to a photometric catalog derived by ourselves 
from the GOODS-South HST-ACS data (GOODS data 
hereafter), using SExtractor 3 (Bertin & Arnouts 1996). 
The resulting sample contains 505 objects. 

To analyse the surface brightness profiles of our galax- 
ies in a similar rest-frame band within the explored 
rcdshift range (0.1<z<l.l), we have extracted the pro- 
files in the following bands: Veo6 band for galaxies with 
0.1<z<0.5, 2775 band for 0.5<z<0.8, and z 850 band for 
0.8<z<l.l. This allows us to explore the surface bright- 
ness distribution in a wavelength close to the B-band 
rest-frame. 

3. RADIAL PROFILE ANALYSIS 

The surface brightness (yU-r) profiles were extracted 
using photometry on quasi-isophotal elliptical apertures. 
The intensities were estimated as the median of the flux 
in the area between elliptical apertures of increasing semi 
major axis length (hereafter we will dub these lengths as 
"radii" ) . The ellipticity and position angle of the aper- 
tures were fixed to those retrieved by SExtractor for the 
whole distribution of light of the object. The center of 
the apertures was also fixed: in a first iteration, the first 
moments (in l x J and 'y') of the brightness distribution 

3 This catalog was obtained by detecting sources in the 2850 band 
which had at least 16 contiguous pixels (0.014 arcsec ) at an 
isophotal level of 0.6 sky a (25.35 mag/arcsec 2 ) or higher. These 
photometry parameters are the same employed by the GOODS 
Teams for producing their catalogs. Nonetheless, we wanted to 
have segmentation maps for each object, as produced by SExtrac- 
tor, but the GOODS team did not release these files, and so we 
had to do our own, though equivalent, photometry 



of flux of the object were used for this purpose. After 
visual inspection of the resulting profile and the image 
of the object, the center was refined, if needed, with help 
of the task "imexam" from iraf 4 to match what was vi- 
sually estimated, in each the dynamical center 
of the object. In a regular disk galaxy, as are those un- 
der study, this center coincides with the central bulge or 
nucleus. The radii of the annular apertures was linearly 
increased at constant steps of 1 pixel (0.03"), up to a ra- 
dius which is a 50% larger than the radius of a circle with 
the same area as the isophotal area of the object in the 
0850 band, as given by SExtractor. The isophotal area 
is that covered by the set of connected pixels with in- 
tensities above the detection threshold which constitute 
a detection (in our case 25.4 mag/arcsec 2 in zsso)- For 
objects with a more or less regular morphology, as are 
those selected, the median intensities in these annulii arc 
a good approximation to isophotal intensities. The error 
in the intensity, SI, is given by the a of the distribution 
of fluxes inside the annulus, divided by the square root 
of the corresponding pixel area. The intensities ("/") 
were transformed to surface brightnesses (p), expressed 
in mag/arcsec 2 in the AB system, using the magnitude 
zero points posted on the GOODS-HST/ACS website for 
each filter and the angular scale per pixel, through the 
expression [i — zero — 2.5 • log(I/ scale 2 ). For the er- 
rors in a, Sfi, the formula used was S/i = 2.5 • log( 1 + 
SI /I). In producing these profiles, SExtractor, DS9 and 
Iraf software packages were used, "glued" together by a 
script written by ourselves in "Python" 6 language. 

We produced radial profiles, and characterized them 
in terms of the properties of the exponential laws which 
best fit them. For each object, the image and profile 
obtained are visually inspected. If the object is affected 
by artifacts (which mostly occurs near the edges of the 
area of the GOODS-HST-ACS field), or the object is sus- 
pected of being a merger, or has an irregular morphology, 
or the object is between two tiles of the survey, or part of 
the object lies outside the field, the object is rejected for 
further analysis. From the initial 505 objects, 70 were 
rejected for these reasons (14%). 

For those objects which are suitable for analysis, the 
profile was tested for the existence of breaks of any kind, 
trying the 3 methods we describe next and comparing 
results, and also by visually judging the apparent mor- 
phology of the galaxy. If no break is apparent (Type I), 
the profile is fitted to a single exponential. If there is a 
break, it is characterized by using one of the two follow- 
ing methods. The first, which we call the "Intersection 
method", was the most commonly applied. It was used 
in 84% of the 289 objects which show a "break" (i.e. 
Types II and III). In this method two exponentials were 
independently fitted to non overlapping sections of the 
profile. The point of intersection of the two lines gives 
the Break Radius, Rb t , while the surface brightness on 
the profile, linearly interpolated at that radius, is the 
surface brightness at the break, [i s r ■ The radial scales 
of the exponentials (hi and h 2 ) are used to distinguish 
between Types II (hi > h 2 ) and III (hi < h 2 ). 

In 47 cases (16% of the objects which show a break), 

4 http://iraf.noao.edu 

5 http://hea-www.harvard.edu/RD/ds9/ 

6 http://www.python.org 
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a slightly more involved process is required to measure 
the position and brightness of the break. We have found 
some objects which show a change in slope of some Type 
(mostly II, but also III), in which there is a smaller zone 
between the 2 exponential regions with a different slope 
from the other two more extended zones. An example 
of this phenomenology is shown in Fig. 1 . If the length 
of this intermediate zone were increased, and if its slope 
were constant, it would be classified as a "mixed" Type 
(II+III or III+II, Erwin et al. 2005, PT06), but' we 
refer here to those cases in which the radial extent of 
this intermediate zone is so small that it seems to be 
showing a specific morphology. Measuring the position 
of the break in these (from now on in this paper) 
termed "Jump Truncations", cannot be done with the 
"Intersection method". In these cases, the intersection 
point does not match the point at which there is the 
change of slope in the profile, as seen in Fig. 1. Instead, 
we use the "Equal Deviation" (ED) method, which we 
have devised for this purpose. It consists in locating the 
position of the break at the point in which the inner and 
outer fit lines, are at a maximum and equal distance (in 
/z) from the intensity profile. In most cases, there are 
usually two points of the profile which are at a same 
distance from the fits. One is the point of intersection of 
the inner and outer fit lines (a meaningless position in 
these cases), and the other is the point of interest, the 
point at which the change in slope appears to happen. 
By selecting the point at which the distances are equal 
and of maximum absolute value, the required point is 
usually naturally selected. Only in some cases, with 
more irregular profiles, the region in which we search for 
the points of equal deviation must be limited to find a 
meaningful solution. If we apply this method to a profile 
in which there is normal break, i.e. without a "Jump", 
the same solution is obtained by this method as with 
the standard intersection technique. It is important 
to note that this method is used only to measure the 
radius of the break (understood as the change in the 
slope of the profiles) with greater accuracy. We think 
that, at this point, and for the stated purposes, it is not 
necessary to speculate about the causes for these "Jump 
Truncations" . Our standpoint is to consider them as 
irregularities in the profiles which make the task of 
estimating the position of the "break" more involved 
than in most other cases. In analysing the profiles we 
adopt the following assumption. A "break" in a stellar 
disk is a significant discontinuity in the slope of the 
exponential profile. This significance comes from both 
the change in slope and the permanence of the change 
at increasing radii. Bearing this in mind, the position 
of the break should be placed at the point where the 
change in slope takes place, and not where the lines 
which best fit the two subsections of the profile intersect, 
if these two positions do not coincide, as is the case for 
the "Jump Truncations" . 

Taking into account that stellar disks have their ir- 
regularities and asymmetries, it might be argued that 
perhaps other profile analysis techniques would be bet- 
ter fitted to our goals. For example, using averages of 
profiles in different sectors of the galaxies, instead of the 
complete azimuthal profiles, might be helpful to indicate 
problems in individual extracted profiles. Three are the 



main reasons which make us prefer the adopted method- 
ology First, our profiles are not azimuthally averaged, 
but we extract median values of intensity, so minimizing 
the relevance of asymmetries on the resulting intensity 
profiles. Second, averaging profiles in sectors would cause 
a severe decrease in the limiting surface brightness which 
could be analysed. For many cases, this could imply that 
we would not be able to measure the truncations. And 
third, selecting the sectors could add a bias which could 
compromise the statistical significance of the results. 

The complete list of objects (505), with detailed results 
of the profile classification and characterization, is given 
in Table 1. Also listed are the objects which were rejected 
for analysis. These have no information on classification, 
and the reason for rejection is given in the last column. 
We will discuss the results in the following sections. 

4. COMPLETENESS & STRUCTURAL PARAMETERS 
RELIABILITY 

4.1. Completeness 

The selection criterion employed for the absolute B 
magnitude (Mb < -18.5) is the same as that used in Tru- 
jillo & Pohlen (2005). They adopted this "cut" as a 
compromise between maximizing the number of objects 
in their final sample and assuring as much homogene- 
ity (equally luminous objects) as possible through the 
range of rcdshifts they explored, which is the same as 
ours. In Fig. 2 (left) we see that at redshift z^0.7, 
which is the mean/median redshift of our sample, the 
completeness level is at Mb ~ -19.5 mag, a magnitude 
below the selected level. This makes our sample incom- 
plete at fainter magnitudes beyond z~0.5. Nevertheless, 
we decided to maintain our absolute magnitude selection 
criterion to satisfy the following objectives: a) to allow 
a detailed comparison with TP05, and b) because oth- 
erwise our sample at low redshift would be significantly 
reduced. It is also interesting to note from the same fig- 
ure that our sample lacks some high luminosity objects 
at lower redshifts because the surveyed volume is signif- 
icantly smaller in that range. In the right panel of the 
same figure we see the completeness level in stellar mass 
is M± ~ 3- 10 9 Mq at z^0.7. The incompleteness with 
stellar mass is less severe than with luminosity. This is 
because the M+ /L relation evolves with redshift in this 
range, in the sense that for a same M* galaxies at z^l 
are brighter in the B-band than at z^0 (eg. Rudnick et 
al. 2003). 

4.2. Comparison with Radial Profile analysis in TP05 

To check the accuracy of our structural parameter de- 
termination and galaxy type classification, we made com- 
parisons with deeper observations and with simulations. 
Trujillo & Pohlen 2005 (TP05) studied radial profiles for 
a sample of 36 late-type galaxies, using imaging data 
from the Hubble Ultra Deep Field (UDF; Beckwith et al. 
2006), which is also contained in GOODS-South. Those 
objects are included in our sample, and 35 of them were 
analysed. Only one object (UDF3203) was discarded in 
this work, as it has B = -18.46 > -18.5 mag, which is 
one of our selection criteria. Of the 35 matched objects, 3 
more objects were discarded in our analysis: one because 
it was near the edge of a tile (UDF7556), and the other 
two as being probable merger candidates (UDF8049 and 
UDF8275). In Fig. 3 we show some examples of radial 
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profiles using UDF data and GOODS-S data. We see 
that the profiles match well (\hudf — HGOODs\ ^ 0.2 
mag/arcsec 2 ) out to a level of ~ 26 mag/arcsec 2 , and 
then the differences grow erratic, as the UDF data are 
"deeper" . We take this value as a minimum surface 
brightness level of reliability for our profiles. 

Of the 32 objects which were in the TP05 sam- 
ple that we studied, in 22 cases (69%), we assign 
the same classification for their profile (as Type 1, II 
or III). The remaining 10 cases (31% of the matched 
sample) show varying bases for disagreement. In 4 
cases (13%, UDF3268 IItpos-Wbos, UDF3822 IIItpos- 
I at bos, UDF9455 IItpos-Utbos, UDF6853 Hippos- 
I at bos) the feature which could be taken as a "break" 
can also be identified in our profiles, but happens in TP05 
at a lower surface brightness part of our profiles which for 
our purpose is unreliable (i.e. HBr^ 26 mag/arcsec 2 ). 
In 2 cases (6%, UDF6862 III T po 5 -II ATB08 , UDF 8257 
IIItpos-Hatbos), we have suspicions of neighbor objects 
being the cause for these galaxies being classified as Type 
III in TP05. Three galaxies (9%), present morpholo- 
gies which are rather irregular or have profiles of diffi- 
cult interpretation, because of bars (UDF2525 -II+III?- 
n T P05-IIlATB08, UDF8040 - "Jump Truncation"- I T P05- 
IFatbos, UDF7559 -barred- IIItpos-IFitbos)- And only 
in 1 case (3%, UDF4491 IItpos - Pitbos) we can see no 
feature in our profiles that may justify the difference in 
classification. More importantly for this paper, in the 
case of Type II galaxies there were only 2 cases out of 
19 (^10%) that were misclassified due to an insufficient 
depth of the GOODS images. 

In Fig. 4 we present a comparison of the estimates 
of the i?e r and fi Br ior the 15 objects of Type II where 
our classification agrees with TP05. We see an over- 
all good agreement for both parameters. We define 
the relative error in R Br as 100 * [(-Rtpos-^tbos)/ 
Ratbos], and the relative error in intensity as 100 * [1 
_ iQ-0A(u TPa5 -u AT BO8)^ Now, for the given comparisons, 
the standard deviations are or = 0.2" (crff = 22%) 
and a u = 0.6 mag/arcsec 2 (af r = 43%) for i?e r and 
\iBr respectively. For R Br , the largest disagreement is 
with UDF1971. The feature we identified as the break 
was taken by TP05 as a mere consequence of the bar, 
while the truncation they detect falls below our reliabil- 
ity level, though it can be seen in our profile. In fj, Br the 
differences arc below 1 mag/arcsec 2 in all cases, except 
for UDF6974 and UDF1971. There is a hint, although 
based on poor statistics, that breaks at \i Br (TP05) > 25 
mag/arcsec 2 are overestimated in brightness in our work 
relative to TP05 {/j, Br (ATB08) < [i Br (TP05)). 

4.3. Profile. Analysis of Artificial Galaxies. 

In addition to the previous comparison, we performed 
simulations in order to assess the accuracy of our re- 
sults, with regard to the classification of the profiles and 
their characterization. We have created 300 synthetic 
galaxies using the GALFIT program (Peng et al. 2002, 
see the Appendix for a detailed description of how the 
objects were created and analysed). The mock objects 
were placed in a synthetic image, to emulate the signal- 
to-noise ratio, angular scale and PSF properties of the 
zssoband images of GOODS-S. The synthesized galax- 
ies are exponential disks of Type I, II or III (100 of 



each kind). Their total apparent magnitudes, inner scale 
length h\ over every Type, and outer scale length I12 , 
R Br and \i Br for Types II and III, are taken to cover the 
whole ranges over which those parameters vary for our 
real objects. The selected ranges of parameters are in 
Table 4. 

The simulated objects were detected with SExtractor 
and analysed in an analogous way to that used in study- 
ing the real objects. Each model was visually classified as 
Type I, II or III. When the object had been assigned to 
Type II or III, the position of the break in the /x-r plane 
was also estimated. The process was done "blindly" , the 
classifier did not know which were the real classification 
and profile parameters of the object a priori. 

In Fig. 5 (left) we show the effect of the apparent mag- 
nitude of the object on its classification. In the panels 
we have plotted the fraction of objects of each Type (I, 
II or III) which are correctly classified, against the total 
apparent magnitude. In each panel the lower 2 curves 
represent the fraction of objects which are misclassified 
as being of the complementary Types. Only points which 
represent a population of at least 3 objects are shown. In 
the first panel we see how the success ratio in classifying 
Type I objects goes from 100% to 50% at the fainter end. 
The wrong classifications are divided between Types II 
and III with roughly equal shares. In the second panel we 
see how Type II objects are classified with a higher ratio 
of success at all surveyed magnitudes (>70%), but fainter 
objects are principally misclassified as being of Type I. 
Finally, as shown in the third panel, Type III (>80%) 
objects are never misclassified as Type II, but Type I. 
It is interesting to note that the reliability of our clas- 
sification, both in models and real galaxies (Sect. 4.2), 
are at a similar level: 70-80%. Furthermore, note that 
more than 80% of our real galaxies have z 85n <23 mag, 
for which the probability of a potential misclassification 
is lower. 

As an illustrative exercise, we have also tested the re- 
liability of the visual appreciation of the position of the 
"break" in the simulations. In Fig. 5 right we repre- 
sent the relative error in the estimate of Rs r , against 
the input surface brightness at the "break" fi Br i n . The 
relative error in R Br is defined as 100 • (R Br out - Rbt in) 
I Rbt- in, where R Br in is the actual position of the break 
for the model and R Br out is the visually estimated posi- 
tion. In the Figure we plot the results for Types II and 
III. The points have been fitted by a line, and this gives 
a variation in the relative error from -9.4% to 3.9% be- 
tween 22 and 26 m&g/arcsec 2 in fiBr in- The standard 
deviation in the relative error is er err = 8%. We see how 
the scatter for Type II objects is less (a^ r —5%) than for 
Type III objects (<7^=13%). These errors are both be- 
low the estimate of the error in R Br (22%) derived from 
the comparison of results in this work and TP05 for ob- 
jects common to GOODS and UDF observations, as re- 
ported in subsection 4.2. It is that error, derived from the 
analysis of real galaxies, which must be taken a realistic 
estimate of the error in the estimate of R Br ■ 

Finally, in Fig. 6 we present probability maps of the 
success in classification of simulated objects of Type II 
(left) and III (right), as a function of the simulated h\ and 
fJ-Br (up) and i?e r and fx Br (down) of the models. The 
success ratio has been coded in gray scale, and the por- 
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tions of the planes (hi - /J-Br and Rs r - M-Br ) f° r which no 
input model within the corresponding ranges of param- 
eters was produced are marked with a cross. The points 
mark the position in the given planes of the real galaxies 
under study. Most of our real galaxies are in a position 
of the plane where the classification success is >80%. 

To summarize the results of this section, from a com- 
parison with deeper observations (i.e. with the UDF) 
and with artificial galaxies in simulations we conclude 
that for Type II galaxies (the main goal of this paper) 
the accuracy on type profile identification is higher than 
80%, for surveys as deep as GOODS-South. Nearly 20% 
of Type II galaxies that are potentially missed in our 
work are not found because their break positions lie at 
fainter magnitudes than 25-26 mag/ arcsec 2 . 

5. RESULTS 

5.1. Classification of Galaxies in our Sample. 

We have divided our sample in 3 redshift bins or ranges; 
"low": 0.K z <0.5, "mid": 0.5< z <0.8 and "high": 
0.8< z <1.1. Of the 435 objects which were suitable for 
analysis, 146 (34%) are classified as Type I, 242 (56%) 
as Type II and 47 (11%) as Type III. All considered ob- 
jects are visually confirmed as non Irregular/Merger, as 
we stated before. There are 6 objects, of Types I and 
III, which are suspected of being of early-type (probably 
SO, though it is difficult to be sure without more infor- 
mation), but have entered our statistics. These objects 
are marked as "Early" in Table 1, and the reason to in- 
clude them is that they seem to have a genuine disk. 
Nonetheless, and to avoid any controversy on this issue, 
our results on disk sizes are based on a subsample of the 
Type II objects, of which none is classified as "Early". 

Our work is focused on the population of late-type, non 
interacting galaxies, in which the profile is "externally" 
truncated. This means we concentrate on Type II ob- 
jects, and specifically on those in which the "downbend- 
ing" break takes place in the outer parts of the visible 
disk, i.e. beyond the visible bars or spiral arms. A se- 
lection criterion for the ratio Rbt / hi , where hi is the 
scale length of the inner part of the disk, is not effective 
in discriminating between the two kinds of "breaks" : the 
"inner" , which take place "inside" the stellar disk, and 
the proper "truncations" , those which mark the edge of 
the stellar disk, or a decline in the density of the stellar 
population, and so the trained judgement of the classifier 
is needed to distinguish between them. We define, then, 
our "truncated" sample (T-sample hereafter) as follows: 
late-type objects, though not Irregular, nor involved in 
merger events, with a Type Il-truncated profile. This 
sample contains 238 objects. That means only 4 out of 
242 objects which form the whole Type II sample have 
breaks of the "inner" kind. Additionally, 48 of the ob- 
jects in the T-sample have asymmetries in their disks 
(20%). Finally, there are a few objects (15, 6.3%) in 
which the profile is of the mixed kind II+III. In these 
there is a clear Type II break, followed by a second break, 
of Type III; even so we have included them in the T- 
sample. 

The distribution of the objects amongst the 3 redshift 
ranges, is given in Table 5. We find that the frequency 
of objects classifed as of Type I decreases from 39% to 
25% between z~l and z~0.3. For Type II objects the 
opposite trend is found, increasing their proportion by 



almost 9% in the same period. These results are com- 
patible with no changes in the distributions within the 
error bars. Finally, Type III objects show no change in 
their relative population within the error bars between 
z^l and z~0.3. It is tempting to explore whether there is 
a real evolution in the populations of objects of different 
profile types. This would mean that more disk galax- 
ies presented no truncations in the past than nowadays, 
relative to the fraction of those that present truncation. 
But first we need to look again Fig. 5 (left), in which we 
show results from the simulations. As we stated above, 
when Type I objects are fainter, they are increasingly 
mistaken for Types II or III in roughly equal proportions. 
On the other hand, genuine Type II objects are increas- 
ingly classified as Type I when the magnitude decreases. 
Consequently, it is quite possible that the observed vari- 
ation in the frequencies is caused by misidentifications of 
Type II galaxies as Type I at high z. 

We also compare the given ratios of profile types with 
results published in PT06 for local disk galaxies in the 
SDSS (their sample being composed of 85 objects, with 
Hubble Types T in the range 3<T<8.5). In the Table 4 
of that paper they report the following shares: 11±3% 
of objects are classified as Type I, 66±5% objects have a 
Type II profile and 33±5% have a Type III profile. The 
shares do not add up to 100, and that is because they 
count objects with a mixed classification (i.e. II+III or 
III+II) for the statistics of both Types. In our lowest 
redshift bin (z^0.3) the objects are classified as 25±6% of 
Type I, 59±9% of Type II and 15±5% of Type III. We see 
the fraction of Type II objects is in good agreement with 
their results, while there is significant disagreement with 
regard to Types I and III. We assume that the fact that 
Type III is under represented in our work, in comparison 
with theirs, is probably because Type III objects have 
their breaks at lower brightness levels than Type IPs, 
and so are harder to identify at intermediate redshifts. 

5.2. Results for the "truncated" sample. 

As stated in Section 2, the classification and charac- 
terization of the profiles has been performed in the band 
which best approximates the rest-frame i?-band in each 
redshift bin. As a consequence, the contribution of lu- 
minous young stars to the retrieved surface brightness 
profiles is significant. In fact, in many objects the sig- 
nature of star forming "clumps" spread over the disks 
is evident in the images. This means that the reported 
"truncations" should be more accurately interpreted as 
related to, or at least influenced by, the extent of the 
"star formation" disk, a term we will use to refer to that 
part of the disk where most of the massive star formation 
is taking place in the galaxies. In contrast, the "trunca- 
tions" found at longer wavelengths in other works (e.g. 
PT06) could be more unambiguously interpreted as an 
abrupt decrease in the density of stellar mass. Nonethe- 
less, when comparing profile parameters of objects in the 
Local Universe (from results in PT06) and at intermedi- 
ate redshift from this work, we do so based on profile 
characterizations performed in similar rest-frame bands: 
g' for local objects and i?-band for the higher redshift 
ones. 

5.2.1. Rbt - B- Luminosity relation. 
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Using the radial position of the truncation as a direct 
estimator of the size of the "star formation dominated 
disk", we have explored the relation between i?s r and 
the B-band Luminosity of the galaxies for the T-sample 
("truncated" profiles). In Fig. 7 left we show Rsr in kpc 
against the rest-frame absolute B magnitude (from Bar- 
den et al. 2005). In the 4 panels we show this relation at 
the 3 redshift ranges we have explored, and also in the 
Local Universe, based on results from PT06 on SDSS 
galaxies at z^O (we use the values of Rs r obtained in g'- 
band for these objects). The observed distributions have 
been fitted to a line with a "robust" least absolute devi- 
ation fit, in an iterative "bootstrap method" to get the 
most probable values for the slope and the y-intercept 
ordinate and their errors (the same mathematical pro- 
cedure is used in Figs. 8 and 11 through 14). At a 
given luminosity, Rs r is evolving towards smaller values 
as redshift increases. 

Before exploring the meaning of this last assertion we 
want to clarify an important point with regard to the 
distribution of objects with luminosity at different red- 
shifts. In the panels we see how this distribution varies 
amongst redshift bins. In the "low" redshift bin the ob- 
jects are concentrated around (-20.0, -18.5) mag, while 
in the "mid" range the distribution is quite homogeneous 
between -18.5 < Mg < -22 mag. Quite opposite to what 
happens at lower redshifts, in the "high" redshift bin the 
objects cluster around (-20.0, -21.5) mag. These distri- 
butions are the result of a) probing a smaller volume in 
the lowest redshift bin and b) the sample is incomplete 
at the lower end of luminosities in the higher redshift 
bin. To avoid as much as possible the sampling effects 
just described, we explore the size evolution using as ref- 
erence point M B = -20 mag, which is well populated of 
galaxies at all redshifts. 

In Fig. 7 b) we show the value of the aforementioned 
best fit lines at Mb — -20 mag, relative to the corre- 
sponding value in the z—0 sample (logarithmic y-axis), 
against the mean value of redshift in each bin. The er- 
rors shown are from the bootstrap method applied to the 
fitting of the relation. We see how the ratio decreases 
with redshift. Given the fact the points are aligned we 
fit them with a straight line. On this line, the Rbt of 
galaxy with Mb = -20 mag has increased by a factor 
2.6±0.3 between z=l and z=0. We see how the point at 
z=0.3 fits somewhat worse to the straight line than the 
others. Nonetheless, the deviation is 2.5o\ 

As we explained above, the ranges of Mb covered by 
the objects are not the same amongst different redshift 
bins, because of a lack of completeness inherent to the 
available data. Although our analysis method has been 
devised to minimize the negative consequences this fact 
could have on the produced results, it would be desir- 
able to test, within the possibilities granted by the avail- 
able data, whether this method is really effective for the 
stated purpose. In this line, we have also explored the re- 
lation between Rs r and Mb when a same range of Mb is 
used to select the objects in all redshift bins. This range 
has been chosen as -19. 5> Mb >-21 mag. When apply- 
ing this restriction, the sample populations reduce to 15, 
22, 75 and 41 objects in the z~0, and the "low", "mid" 
and "high" redshift bins respectively. In comparison, the 
whole samples explored in Fig. 7 (and also in Figs. 8 
through 14) are 39, 39, 133 and 66 objects in the same 



redshift bins, i.e., significantly larger. The growth factor 
in Rbt between z~l and z^0, deduced from the best fit 
lines to the relation Rs r -Mb at Mb =-20 mag is, when 
restricting the range in Mb , 2.0±0.8, which is smaller 
than the value given above, but in agreement within the 
error bar. This further supports the hypothesis that the 
reported change in Rsr at a fixed Mb across redshift is 
not due to differences in the luminosities range of the 
sampled objects. 

5.2.2. Rbt - Stellar Mass relation. 

It is not entirely straightforward to interpret the size 
evolution of galaxies from the Rs r - Mb relation, studied 
above, for the following reason. The Mass Luminosity re- 
lation, specially for shorter wavelength bands such as the 
B-band, has varied significantly between z=0 and z~l 
(Brinchmann & Ellis 2000; Bell ct al. 2003; Dickin- 
son et al. 2003; Perez-Gonzalez et al. 2008). Objects 
of the same mass were brighter in the i?-band in the 
past as their stellar populations were younger on aver- 
age. This has as consequence that even if there were no 
changes in the i?s r of objects, their proved increase in 
£>-band Luminosity would make the relation vary such 
that objects of a same Mb would have smaller values of 
Rbt in the past. This is the trend reported in the preced- 
ing subsection. Fortunately there is a way to overcome 
the consequences of this luminosity evolution in order to 
test for a "real" change in the sizes of stellar disks with 
time (or more properly, the size of the "star formation 
dominated" disk, as we are relying on rest-frame -B-band 
data). This is to survey the Rbt- M± relation. It is 
also known that the stellar mass content of galaxies has 
increased since z~l, but this change has been more mod- 
erate, in relative terms, (Rudnick et al. 2003, <30%; ) 
than that in luminosity. 

Following the preceding discussion, we have tested the 
relation between size (Rbt- ) and stellar mass. This is the 
strategy adopted in Trujillo ct al. (2004) and Trujillo 
et al. (2006) to test for evolution of the effective radius 
of massive galaxies (M* > 10 10 M Q ) at high redshifts. 
In Fig. 8 we show, in analogous fashion to Fig. 7 left, 
the values of Rbt for objects of the T-sample against the 
stellar mass M+ , as reported in Barden et al. (2005). 
The z=0 data come again from Trujillo & Pohlen (2005) 
(their stellar masses were computed following the pre- 
scription of Bell et al. 2003). In this case, the best 
fit line (obtained by the same method as in the Rbt~ 
Mb relation) also does show a trend to lower values with 
redshift. While in this case the differences in the distri- 
bution of stellar masses are not as large as in the case 
of the luminosities, we see how the slopes are slightly 
less "stable" than in that case, based on the assumption 
that they should follow less pronounced changes between 
redshift ranges. This may be because the distribution 
of masses is not as broad relative to the distribution of 
Rbt as the distribution of Mb 's. Moreover, the stellar 
mass is a quantity derived from the luminosities, and 
this leads to an increase in the dispersion. The factors 
that go into producing the dispersion in the two relations 
[Rbt - Mb and i?s r - M*), are the following: a) inaccu- 
rate redshifts (5z/(l+z)~0.02), b) inaccurate estimates 
of the Rsr (cfj" r iS 25%), and c) intrinsic dispersion in 
the relation. Though the two first factors (especially the 
second) may play a role in explaining part of the disper- 
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sion, we do not rule out the third possibility; i.e. it is not 
only possible, but almost to be expected, that objects of 
the same mass may have followed different evolutionary 
paths (depending on initial conditions and environment, 
for example), which have shaped their structure in a dif- 
ferent way, and so produce "intrinsic" dispersions in the 
Rbt- Mb and Rbt~ relations. 

In Fig. 8 right we present the fitted value of Rs r at 
10 10 M Q relative to the z~0 value, against the mean 
value of redshift in each bin. The errors are again de- 
rived from the bootstrap method applied when fitting 
the Rbt -M+ relation. We see a less strong evolution with 
redshift than for the Rsr -Mb relation. It is also inter- 
esting to see how in this case the points are much better 
aligned than for that relation. 

From the linear fit to the points in Fig. 8 right we de- 
duce an increase by a factor of 1.3±0.1 in Rs r for a given 
stellar mass between z=l and z=0, i.e. in the last ~8 
Gyr using the standard parameters in our cosmological 
model. For comparison, in the same range of redshifts, 
and using the assumption of a maximum evolution in 
luminosity for their galaxies, Trujillo & Pohlcn (2005) 
reported a more moderate growth in i?s r by a factor of 
at least 1.25, which is in good agreement with our result. 

Another point worth noting is that individual galaxies 
also evolve in stellar mass. For this reason it is interest- 
ing to provide a rough number of how much Rs r could 
grow for a given "individual" galaxy. In the range of ages 
probed here it is claimed in the literature that galaxies 
have increased their stellar masses by ^30% (Rudnick et 
al. 2003). This means, that a Type II object with M* = 
10 10 Mq at z=l, would have its R Br larger by <50% to- 
day, i.e. by somewhat more than if its mass were un- 
changed. 

We have also explored the relation between Rs r and 
M*when a same range of M*is imposed to select the 
objects in all redshift bins, in analogous way to what 
we did in subsection 5.2.1 for the Rs r -Mb relation, and 
for the same reasons. In this case, the chosen range of 
stellar masses is 5 • 10 9 < M± < 5 • 10 10 M© . In this case 
the populations reduced to 30, 11, 59 and 40 galaxies 
within the z^0, and "low", "mid" and "high" redshift 
bins, respectively. As a result, and for these stellar mass- 
restricted samples, the growth factor in Rs r between z~l 
and z~0 is 1.4±0.2, a minor difference with the value 
obtained for the unrestricted samples. This difference 
between "restricted" and "unrestricted" values is lesser 
than that found for the Rs r -Mb relation because the 
distributions of M* are more similar between the redsfhift 
bins than the distributions in Mb ■ Finally, this test 
strengthens the significance of the result found for the 
growth in the Rbt -M* relation, as it persists when only 
objects within a same stellar mass range are taken into 
account through redshift. 

5.2.3. Surface Brightness at the Break Evolution. 

We have also explored the distribution of /j Br in our 
sample of truncated profiles. In Fig. 9 (left) we show 
a histogram of the fiBr distributions of the T-sample in 
the 3 redshift bins under study. The /i Br have been cor- 
rected for the cosmological dimming effect (I oc (l+z)~ 4 ). 
This means we are showing the distribution of fiBr as 
they would be measured from a rest-frame observational 
standpoint for every object. Also in the same panel 



we represent the median of the distribution of \xb t in 
the g' band reported in PT06 for z^0 galaxies. We 
use this band as it is the closest to our rest-frame B- 
band. We see a clear evolution in the median values of 
the distributions, in the sense that the "break" in the 
profiles happens at a brightness level which is 3.3±0.2 
mag/arcsec 2 brighter at z~l than at z^0. This is a 
strong evolution, by a factor of 20.9±4.2 in intensity. As 
we are measuring the surface brightness profiles in the 
rest-frame £?-band, which is significantly affected by the 
contribution of young stars, we argue this change may 
be related to the well known cosmological evolution in 
global SFR between z~l and z^0. Medians and stan- 
dard deviations for the distributions represented in Fig. 
9 are given in Table 2. 

If we apply here restrictions both in Mb (-19. 5> 
M B >-21 mag) and M* (5 • 10 9 < M* < 5 • 10 10 M Q ) to 
select the objects in all redshift bins the results do not 
change perceptively either, as in previous cases. The dif- 
ference in median values of /iBr between z~l and z^0 is 
-3.5±0.3, only slightly larger than for unrestricted sam- 
ples. This corresponds to a decrease in intensity by a 
factor 25.1±8.0. In this case the samples reduce to 11, 
11, 39, and 19 galaxies in the reshift bins termed as "lo- 
cal", "low", "mid" and "high". 

5.2.4. Rbt / hi Evolution. 

We also present results on i?s r , relative to the scale 
length of the inner exponential, hi , in our truncated ob- 
jects (T-sample). In Fig. 10 are shown histograms of 
the distribution of this parameter in the 3 explored red- 
shift bins. Also represented is the median value of the 
ratio for the local sample in PT06, Rbt /hi (z~0) = 2.0 
(measured in (/-band) . The statistical parameters of the 
distributions can be found in Table 3. The most striking 
feature is the low set of values of the ratio Rs r /hi we 
have measured, compared to local values. The mean, 
median and standard deviation of the distributions of 
values in the different redshift bins are given in Table 
3. The median values are also shown as vertical lines in 
Fig. 10. Perez (2004) found a median/mean value of 1.8 
for Rbt /hi at z~l (over 6 objects), which is also shorter 
than the local value, but still larger than those we find. It 
is probable that the difference in the sizes of the samples 
is the cause of this difference. With regard to the dif- 
ference between the values of RBr /hi she obtained and 
those observed in the Local Universe, she attributed it 
mainly to two biases: 1) the detection limit on surface 
brightness would have prevented her from identifying the 
breaks with larger value of Rs r /hi ; and 2) the effects of 
dust. The first reason is evident, but the latter may need 
clarification. She assumes dust absorption was more im- 
portant in late- type galaxies at intermediate redshifts, 
and specially in the inner parts of the galaxies. This 
would be perceived as the inner parts of disks having 
their profiles "flattened" , i.e. with larger hi . This is 
what would make the quotient Rs r /hi smaller at higher 
redshift, relative to the Local Universe. Again in Fig. 
10, it is also interesting how the distributions of the ra- 
tio Rbt /hi&ve quite similar amongst the sub-samples 
between z~0.3 and z~l. 

It is important to note that the Rsr /hi evolution pre- 
sented here is measured in relation to the observed (un- 
corrected for dust) Rsr/hioi the local galaxies; con- 
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sequently, if the dust opacity were not to change with 
rcdshift, the observed evolution presented here would re- 
flect an intrinsic evolution in the profile morphology of 
the objects. However, it is likely that the opacity of the 
galaxies changes with redshift. At a fixed inclination, 
bulge-to-total flux ratio, and rest-frame wavelength, the 
degree of attenuation and the increase in the observed 
scale length due to dust can be parameterized by the 
change in the central face-on optical depth. The optical 
depth is a very uncertain quantity (even in the nearby 
Universe), and this makes a detailed evaluation of the 
effect of dust to go beyond the scope of this work. Con- 
sequently, we have not made any attempt to correct our 
results for the effect of opacity. Nevertheless, in order to 
provide a crude estimate of how a significant increase in 
opacity could affect our results, we have performed the 
following exercise: let us assume a mean inclination of 30 
degrees and an increase in the total central face-on opti- 
cal depth in the B-band from 4 (present-day galaxies) to 
8 (high-z galaxies). This change implies a transition from 
an intermediate to a moderately optically thick case. In 
this case, for a disk-like galaxy observed in the £?-band 
rest-frame, the attenuation increases by 0.25 mag (Tuffs 
et al. 2004, , their Fig. 3 and Table 4) and the scale 
length increases by 12% (Mollenhoff et al. 2006). If 
we account for these numbers, the galaxies in our high-z 
sample would be intrinsically brighter by 25% and their 
scale length intrinsically smaller by 12%. In this sense, 
if we assume that the break position is not affected by 
the dust content, the observed (uncorrected for dust in- 
crease) Rsr /hi evolution presented in this paper would 
be an upper limit to the actual evolution. It is impor- 
tant to stress however, that a ~15% bias in Rs r /hi due 
to dust is far from explaining the amount of evolution 
observed here (Rsr /hi increasing by a ~50% between 
z~l and z~0, as shown in Fig. 10 and Table 3). Finally, 
there is no clear evidence that the mean optical depth in 
dust was in fact so much higher in previous epochs, and 
if the opacity were smaller in the past, then the situation 
would be reversed, with our current estimate of the size 
evolution being a lower limit. 

It could be argued that the low values of Rs r /hi found 
at intermediate redshifts could be related to the pres- 
ence of stellar bars, as objects which present them tend 
to have lower values of this ratio. In fact, in Pohlen 
& Trujillo (2006) it was reported that in galaxies with 
bars (Type II.o-OLR) Rs r /hi is, as measured in r'-band, 
around 1.7 compared to 2.5 in galaxies without bars 
Type II-CT. However, using the largest up-to-date sam- 
ple of face-on galaxies (0.2<z<0.84) from the COSMOS 
2-square degree field, Sheth et al. (2008) find that the 
fraction of barred spiral galaxies declines with rcdshift. 
For galaxies with M* >1O 1O M0 the fraction of bars drops 
from 65% in the local universe to 20% at z~84. It must 
also be taken into account that the Type II.o-OLR ob- 
jects are included in our "local" sample, as, due to insuffi- 
cient angular resolution, it is not possible to discriminate 
with the same accuracy in the higher redshift samples, as 
it is possible to attain with the local objects, whether the 
"truncations" are of this Type (II.o-OLR) or not. These 
facts together imply that the systematical decline in the 
Rbt /hi ratio we see with redshift is unlikely to be due 
to a bar fraction increase. 

Another point which is interesting to discuss with re- 



gard to this problem is related to the apparent similarity 
of many of the high-z galaxies under study in this work 
with the irregular galaxies in the Local Universe. More- 
over, the values of Rs r /hi and hi are similar to those 
derived in Hunter & Elmegreen (2006) for a sample 
of irregular galaxies. However, it is unlikely that the 
present-day irregular galaxies could be the final stage of 
the galaxies in our sample. Our high-z galaxies are much 
more massive than these local irregular galaxies. An in- 
teresting feature, however, shared by these two families 
of galaxies is a large value of their specific star forma- 
tion rate sSFR (see for example the comparison between 
the sSFR fo 10 8 M© galaxies at low redshift and the 10 10 
Mq galaxies at z~0.85 in Fig. 2 of Bauer et al. 2005) 
and the fact that their star formation is well spread over 
the galaxy disk. This is probably the main reason for 
their large inner scale-length. 

With the purpose of shedding some light on this 
puzzling phenomenon, we have also probed how 
Rsr /hi relates to global properties of the galaxies. In 
Fig. 11 (left) we present this ratio against Ms for the 
T-sample, in analogous way to Fig. 7. Again, we have 
divided the objects in 3 redshift bins, and the local data 
are from PT06 ((/-band results). In the low rcdshift bin 
(0.1< z <0.5) our derived relation depends less on lumi- 
nosity than that reported by PT06. In the "mid" red- 
shift bin, though, both relations show as parallel. In 
the "high" range the relation reverses, in the sense that 
Rsr /hi decreases slightly with the i?-band luminosity of 
the galaxies, in contrast to the sense in the Local Uni- 
verse. This feature has not been previously reported and 
we must try to explain it. First, in the "high" redshift 
bin there is a lack of low luminosity galaxies which might 
make the slope of the best fit line tilt downwards. But, if 
this were a real phenomenon (and we will give more ev- 
idence below supporting this), and not an observational 
effect, what would it mean? Assuming an exponential 
profile, the ratio Rs r / hi is proportional to the difference 
between the surface brightness of the inner exponential 
measured at the break and at the center: i?s r //iioc ( 
H-Br - M r =0) )• So what we see is that, at z^l the more 
luminous objects had a profile which changed less in sur- 
face brightness, in absolute terms, from center to the 
break, than the less luminous objects (i.e. the surface 
brightness profile is flatter for the more luminous galax- 
ies). Flat profiles, in these rest-frame wavelengths, imply 
that the star formation was taking place more homoge- 
neously across the whole galaxy disk (a good example is 
UDF3372 shown in Fig. 3). 

In the right panel of Fig. 11 we can see the ratio 
Rbt /hi relative to the local value (PT06), at Mb=- 
20 mag, for every redshift bin, taken from the linear 
fits given in the left panel. We see how there is lit- 
tle evolution in the ratio from z^0.3 to z~l. The ra- 
tio R Br /hi varies between 0.62±0.04 and 0.66±0.11 of 
the local value at these redshifts, although the error bars 
are large to derive any firm conclusion. The points in 
the panel give the impression that there is a "gap" be- 
tween the values at z~0 and those at intermediate red- 
shifts, perhaps raising concerns about some kind of differ- 
ence between samples, other than genuine evolutionary 
changes in this structural parameter, or variations in the 
methods employed to study the samples. We must say 
first that, if we take into account the error bars, the sig- 
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nificance of this "gap" is significantly reduced. We also 
want to stress that the analysis methodology employed 
here is analogous to that used in studying the local ob- 
jects. But we will come over this point below, where the 
relation Rs r -M+ is broached. 

We have also tested how the Rbt /hi for the T-sample 
relates to the stellar mass of the galaxies, and this we 
can see in Fig. 12. In the left panel, we see the given ra- 
tio against the stellar mass for the different redshift bins 
(local data are from PT06, as obtained in the (/'-band). 
Both panels are analogous to those in Fig. 11, but sub- 
stituting Mb for M* . In this case, the relation between 
Rbt /hi and in the "low" redshift bin is steeper than 
the relation reported by PT06 at z~0. Note that the lo- 
cal sample covers a smaller range in stellar mass, which 
could somehow bias the local Rs r /hi determination. In 
the mid redshift bin, the two relations are almost par- 
allel, and in the high redshift range, the slope changes 
from positive to negative, as happened in the Rs r /hi vs. 
Mb in the same range of redshifts. It is important to note 
that the dispersions are quite large, and so these slopes 
must be taken with caution. But it is interesting to note 
again an anticorrelation between Rs r /hi and M* at z~l. 
As we said above for the more luminous objects, it seems 
that more massive galaxies at high redshift have an inner 
disk whose surface brightness decreases less between the 
central part and the break, compared to the less massive 
ones. 

In the right panel of the Figure 12 we show the values 
of the best fits to the distributions of Rs r / hi at M* = 
10 10 Mq , against redshift. This time the decrease in 
Rbt /hi is progressive, and we find that the ratio has 
increased by a factor 1.6±0.3 between z~l and z—0, from 
the linear fit. 

A point of particular interest is that the points in the 
right panel of Fig. 12 appear better aligned than in the 
right panel of Fig. 11. This could mean that the ap- 
parent discontinuity in the Rbt -Mb relation commented 
above may be, at least in part, caused by the signifi- 
cant evolution in B-band luminosity in the same redshift 
range. But we would not like to push this interpreta- 
tion, as the error bars are large enough to be compatible 
with a more progressive evolution of Rs r /hi at a given 
Mb with redshift. 

Also for the Rs r /hi parameter have we tested its rela- 
tions with Af^and M* when the corresponding ranges 
to select the objects are fixed, as in previous subsec- 
tions. In the case for a fixed Mb (selection range: - 
19. 5> Mb >-21 mag), the ratio Rbt /hi at intermediate 
redshifts varies between 0.50±0.08 and 0.59±0.17 of the 
local value. This is a slightly larger difference with local 
values than for the unrestricted samples. This does not 
affect in a perceptible way our discussion on this param- 
eter. For the relation between Rs r /hi and M* (selection 
range: 5 • 10 9 < M± < 5 • 10 10 M Q ), the result re- 
mains just unchanged with respect to the unrestricted 
case. This is, there is an increase in the ratio also by a 
factor 1.5±0.3 between z~l and z=0, from the linear fit. 

5.2.5. Scale length of the disk inner to the Break: hi . 

In a further attempt towards clarifying why we find 
a distribution of i?s r //iiwith lower mean and median 
values at higher redshift than in the Local Universe, we 
have also probed the relation between hi and Mb and 



M^ for the "truncated" galaxies under study. 

In Fig. 13 we present hi as a function of Mb for galax- 
ies of the T-sample, again for different redshift bins, as in 
previous similar figures (7 and 11). In the left panel we 
have also included the relation between "size" (i.e. scale 
length derived from the effective radius as explained be- 
low) and absolute magnitude found by Shen et al. (2003) 
for local, Sersic n <2.5 galaxies in the SDSS, represented 
by their best fit curve (slash-dotted line). Shen et al. 
(2003) published the relation between the Sersic effective 
radius of the galaxies (i.e. the radius that encloses half 
of the light of the Sersic model that best fits to the dis- 
tribution of light of the object), as measured in r 1 band 
and absolute magnitude (see Fig. 6 in their paper). Our 
profiles are retrieved in approximately rest-frame £?-band 
(~ g' band), but we would not expect dramatic differ- 
ences in the size estimates of disks between the two bands 
(cither from scale lengths or truncation radii; see eg. val- 
ues of these parameters in both bands in Table in PT06). 
Assuming an n=l Sersic profile (exponential profile), we 
converted the effective radii (Rs, e ff) to scale length of 
the exponential profile using the expression h = Rs, e ff 
/ 1.676. As their cut-off for classification as a late-type 
galaxy was as broad as n <2.5 (the same as ours), this 
implies that most of their objects do not obey the ex- 
ponential law, and we should term this h as the equiva- 
lent scale length. We still show their curve for reference. 
First, we see how the values of hi in the PT06 sample fall 
above the Shen et al. (2003) curve (hi values from PT06 
were obtained in the (/-band). This result is expected 
since fitting the whole profile of a truncated galaxy will 
produce smaller values of scale length than for the inner 
part. Our results are divided amongst the other 3 sub- 
panels as a function of redshift, and we see how the best 
fit lines to the distributions have much steeper slopes 
than for the local sample. Our galaxies cover a wider 
range in hi relative to local ones, reaching larger limits, 
especially for the most luminous objects. We see in the 
right panel the hi values from the best fits to the distri- 
butions at Mb =-20 mag, as a function of redshift. It 
can be seen how there is a moderate decrease in hi with 
redshift. 

In Fig. 14 (left) we show the relation between hi and 
the stellar mass for the galaxies of the T-sample, at dif- 
ferent redshifts. The local data, as before, are from PT06 
(results obtained in (/-band) . The slash-dotted line rep- 
resents the relation between equivalent scale length, h 
= Rs, e ff / 1-676, and stellar mass given in Shen et al. 
(2003) (taken from their Fig. 10). That is a best fit to 
the relation they find for late-type (n <2.5) galaxies in 
SDSS, in the z' band. Again, we would not expect dra- 
matic changes in size parameters of disks between B and 
z 1 band, and this is the best suited result for comparison 
they give. It is evident that as in the hi - Mb relation, 
the values of hi fall systematically above the Shen et al. 
(2003) relation. The reason for this is the same as the 
given previously: the inner scale lengths in "truncated" 
galaxies are systematically larger than the scale-length 
of the whole disk. We can see how our distributions of 
hi are slightly above (reaching higher values) those for 
the local sample, and their best fit lines have steeper 
slopes than in the lowest redshift range. In the right 
panel of the same Figure we see the best fit value to the 
hi distributions at M* = 1O 1O M , relative to the local 
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value, against rcdshift. The values of hi at intermedi- 
ate redshifts are larger than the local ones by a factor 
that goes from 1.29±0.12 to 1.20±0.11 between z -0.3 
and z —0.9. This overall slight increase in h\ , com- 
bined with the decrease of i?s r are responsible for the 
reported decrease in Rs r /hi at higher redshift. The de- 
crease by a factor 1.3 in Rs r (at a fixed stellar mass), 
times the increase in hi by —1.25 (also at a given stellar 
mass) gives 1.6, in agreement with the reported decrease 
between z —0 and z — 1. From the lack of evolution of 
hi at intermediate redshifts seen in the right panel of Fig. 
14, and the previous discussions on the Rbt~ M+ and 
Rbt /hi relations, we can also say something else. The 
disks have become fainter (in rest-frame B-band), while 
keeping their slopes (ocl/hi ) roughly constant or slightly 
increasing, and at the same time increasing their Rb t , 
as the Universe has grown older (in the surveyed range 
of redshifts) . 

We have also probed these relations under restric- 
tions in the selected Mb (-19. 5> Mb >-21 mag) and 
M± (5 • 10 9 < < 5 • 10 10 M e ), depending on whether 
hi is related to luminosity or stellar mass respectively. 
On one hand, for the hi - Mb relation, this varies by a 
factor 1.02±0.2 between z~l and z— 0; i.e. no evolution 
in the surveyed range of redshifts is retrieved (in the 
"unrestricted" case this factor is 1.4±0.5). On the other 
hand, for the hi - M+ relation, the result is the same as in 
the case when no restrictions in M*are applied. This is, 
from the linear fit we obtain a decrease in hi by a factor 
0.8±0.1 between z— 1 and z— when the mass is fixed. 
Again we see how using the same ranges of luminosities 
and stellar masses in every redshift bin does not vary 
significantly the results. 

5.2.6. Robustness of the reported evolutionary changes. 

The previous sub-sections have shown that both the 
break position Rs r and the surface brightness at the 
break radius \iBr have evolved with cosmic time. It could 
be argued, however, that part of this evolution could be 
due to selection effects, since as shown in Section 4 very 
faint breaks can be missed if they were present in the 
galaxies analysed. We have checked whether this is the 
case by exploring how our surface brightness limit could 
affect the observed evolution. 

According to Sec. 4, the profiles are reliable down 
to //— 26 mag/arcsec 2 . Being conservative, we esti- 
mate that we should be able to identify breaks with 
high confidence down to /Us r ,Hmit~25.5 mag/arcsec 2 . 
Using this number as the limiting value we could esti- 
mate what would be the rest-frame (i.e. cosmological 
dimming corrected) surface brightness at the break at 
which the number of objects classified as of Type II will 
start to decline due to a selection effect and not as a 
result of a real evolution. These numbers will be 24.4 
mag/arcsec 2 (z —0.3), 23.3 mag/arcsec 2 (z —0.65) and 
22.6 mag/arcsec 2 (z —0.95). As we can see in Fig. 9, 
these numbers are far away (more than a magnitude in 
all cases) from the position of the peak in the observed 
galaxy surface brightness break distribution. This rein- 
forces the idea that most of the observed evolution is real 
and not caused by this effect. 

Another test we have run is to estimate what would be 
the maximum Rs r that could be measured at a given 
absolute magnitude according to the limiting surface 



brightness proposed above. At the highest redshift bin, 
z —1, the "distance modulus" is 43 mag (k correc- 
tion included, derived from the direct comparison of 
^850 magnitudes to Mb ) • Again, we use a surface bright- 
ness limit of 25.5 mag/arcsec 2 . We consider a truncated 
exponential galaxy with a typical scale length of h = 0.5 
arcsecs. With these values, the largest i?B r that could 
be measured would be: Rbt — 11-1 kpc for galaxies with 
Mb =-20 mag and Rbt =18.5 kpc for galaxies with Mb = 
-22 mag. As can be seen these values are again far away 
from the observational data distribution, indicating that 
most (if not all) of the observed evolution is real and not 
caused by our limiting detected surface brightness. 

6. CONCLUSIONS 

In this work we have presented an analysis of surface 
brightness profiles for a sample of 505 galaxies in the 
GOODS-South field, making use of HST-ACS imaging 
from Giavalisco et al. (2004). These galaxies are all 
classified as disk galaxies, based on the selection of those 
objects with Sersic indexes n <2.5 (Sersic profile fitting 
is from Barden et al. 2005). We have also used publicly 
available data on redshift estimates, absolute rest-frame 
B-band magnitudes and stellar masses from Wolf et al. 
(2001, 2003), Barden et al. (2005). We have classi- 
fied the profiles using the presence of "breaks" in the 
exponential profiles of the disks, and the ratio of inner 
and outer scale lengths (hi /hi ) in Types I (no appar- 
ent break), II ( "downbending break" ; hi /hi > 1) and III 
("upbending break" ; /ii//i2<l). This characterization 
has been performed in bands which, within the explored 
range of redshifts, track the rest-frame B-band. 

We have performed simulations on the classification 
and profile characterization of artificial galaxies, analo- 
gous to that performed on real galaxies. With regard 
to their classification, the worst results are for Type I 
objects, 50 % of which we fail to classify correctly at 
the faintest magnitudes (zsso >22.5 mag AB), while for 
Types II and III the success rates are much higher (>70- 
80%). We have also compared our results on classifica- 
tion and Rsr estimate to those presented in Trujillo & 
Pohlen (2005) for a set of 32 galaxies in common to our 
sample and theirs. Their results are based on the HST- 
UDF observations, which are deeper than the GOODS 
images. Both sets of classifications match in —70% of 
cases. Particularly, in the case of Type II galaxies there 
were only 2 cases out of 19 (—10%) that were classified 
differently in our work due to insufficient depth in the im- 
ages. When we compare the Rs r and \±Br estimated in 
the two studies for Type II objects in common, we find 
no clear bias in any sense. The dispersions are af[ r =22% 
and <7^ rr = 43%, which can be used as estimates of the 
errors in these structural parameters. 

In this work we put special emphasis on the study of 
"truncated" galaxies, i.e. those with Type II profile, 
with the "break" taking place in the outer part of the 
disk. This subsample, which we call the T-sample, is 
composed of 238 objects, an order of magnitude larger 
than any sample in previously published work at inter- 
mediate redshifts. We have studied the relation between 
the radius at which the "break" takes place, Rs r , and 
absolute B-band magnitude and total stellar mass for 
these objects. We find a clear evolution in Rs r with red- 
shift, galaxies with the same values of luminosity/stellar 
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mass having shorter Rs r 's than local galaxies (using as 
reference the results in Pohlen & Trujillo 2006). Wc 
measure an increase by a factor 1.3±0.1 in Rs r between 
z ~1 and z ~0 at fixed M* = 10 10 M Q , and by a factor 
2.6±0.3 at fixed Mb = -20 mag, in the same range of 
rcdshifts. 

At the same time, there is also clear evidence for a 
decrease in the surface brightness level /j Br at which the 
"break" takes place on the profile in the last ~8Gyr. Wc 
find that fiBr (hi the rest-frame -B-band) has decreased 
by 3.3±0.2 mag between z ~1 and nowadays. This is 
equivalent to a decrease in surface intensity by a factor 
20.9±4.2. 

Another point of interest is how Rs r relates to the scale 
length hi of the disk inside the "break" . We show results 
on the ratio Rbt /hi at different redshift ranges up to 
z~l. The median values of this ratio are significantly 
smaller than those found in the Local Universe (1.4±0.6 
at z ~1 versus 2.0±0.7 at z ~0). We also find that 
in the highest redshift bin (0.8< z <1.1) Rbt /hi tends 
to decrease moderately with Mb and M± , in contrast to 
what happens in the Local Universe and the other red- 
shift bins. This means that, at the epoch corresponding 
to z ~1, more massive objects had somewhat smaller dif- 
ference in brightness between the central part of the disk 
and the break radius, than less luminous / massive ob- 
jects (implying that their surface brightness profiles are 
flatter). 

Finally, we have also probed the relation between 
h\ and B-band luminosity and stellar mass. We see how 
the distribution of values of hi is broader, reaching higher 
values, at intermediate redshifts than in the Local Uni- 
verse. The mean values of hi seem to be somewhat larger 
than local ones, by a factor 1.3 at most, but there is al- 
most no evolution in them (within our statistics) between 
z ^0.3 and z ~1. 

We have also tested the previous relations when 
the samples are selected within the same ranges of 
Mb and/or M± in every redshift bin. The reported fig- 
ures do not change significantly under these restrictions, 
and this further supports the robustness of the results. 

We conclude by remarking that our results arc con- 



sistent with the following picture of disk evolution. We 
find that in the lapse of time that goes from z ~1 to 
the present epoch (~8 Gyr), the radii of the "break" in 
truncated disks, as measured in ~_B-band, have increased 
noticeably, while the intensities (on the profile) at which 
the break takes place have substantially decreased. We 
interpret the Rbt , in this particular case, as a measure 
of the size of the disk where most of the massive star 
formation is taking place, as the £>-band is more influ- 
enced by the emission from these young stars than redder 
bands. At a given stellar mass, the scale lengths of the 
disk in the part inner to the "break" were on average 
somewhat larger in the past, and have remained more 
or less constant until recently. This phenomenon could 
be related to the spatial distribution of star formation, 
which seems to be rather spread over the disks in the im- 
ages. So disk galaxies had profiles with a flatter bright- 
ness distribution in the inner part of the disk, which has 
grown in extension, while becoming fainter and "steeper" 
over time. This is consistent with at least some versions 
of the inside-out formation scenario for disks. 
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APPENDIX 

TESTING THE RELIABILITY OF THE CHARACTERIZATION OF PROFILES. 

In this Appendix we explain the work done on the profile characterization of synthetic galaxies, relevant to "breaks" . 
We have generated 300 2D models of galaxies using GALFIT (Peng ct al. 2002), which are simplified emulations of 
real galaxies observed in the range of redshifts we explored. The GALFIT software allows the user to build a model 
as the sum of a series of analytic 2D components. For each component the user can choose the 2D geometry of the 
isophotes, as a generalized ellipse (i.e. with "regular" appearance, but which can also adopt a "disky" or "boxy" 
appearance, depending on the "c" parameter), and besides, also the law that governs the radial profile of intensity. 
The user can choose amongst different laws: exponential, Sersic 1/n profiles, gaussians, moffatians, and more. As we 
are interested in the mid to outer parts of disks, we used only exponential behavior, without introducing bulges or 
other features. It is possible to reproduce a "broken" exponential disk as the composition of 2 concentric exponential 
disks, with the same isophote geometry, cancelling each disk in the internal or external region, relative to the radial 
position of the break. Although we produce these "broken disks" as the combination of two exponential components, 
we then define these models as containing a single component; i.e. a disk with a "break" in its profile. 

All the models synthctized have only one component, with a geometry of the isophotes invariant with radius. We 
chose an elliptical geometry with c=0, i.e. a common ellipse (neither "disky" nor "boxy"). The axis ratio q for each 
model was randomly chosen in the range 0.5<q<l, to mimic variations in orientation of the disks. Three kinds of 
models, each corresponding to a Type of profile (I, II and III) were generated. For Type I galaxies, the model follows 
a single exponential law. The other 2 kinds of models have also a single component, but they follow an exponential 
law with a slope that changes from hi to hi at some radius, the break radius (Rbt), using the technique described 
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above. Depending on the ratio hi j hi , with this kind of model we represent either a Type II (hi / hi > 1) or a Type 
III (hi j hi < 1) galaxy. We generated 300 models, 100 of each of these 3 kinds (Types I, II and III). The ranges of 
h\ , hi, Rbt > HBr and total magnitudes from where to extract the parameters to build each type of model are shown 
in Table 4. The ranges bracket the whole space in which the parameters vary in our sample of real galaxies, at all 
rcdshifts. Some further restrictions are applied to discard certain combinations of parameters, which would produce 
models which are classified as wildly unreal or undetectable by SExtractor: a) the ratio Rbt / hi is limited to be < 6 
; b) the magnitude in the zssoband is within 21<z 850 <24 mag ; and c) the central intensity must be above 3 a s k y of 
the simulated image. 

The 300 created objects were placed in a synthetic, noisy 2; 850 band image with the same pixel scale, zero point and 
noise level as those of the GOODS data. Crowding effects were not simulated as we have discarded in our analysis 
those real galaxies which are too near to other objects. SExtractor was run on the synthetic image, using the same 
detection parameters as those we used for real data, to produce a photometric catalog and segmentation map. These 
were then used to extract the radial profiles of the objects, and we analysed them with the same software we used 
for studying the real objects. Following the same criteria we used for real galaxies, one of us visually classified each 
model as being of Type I, II or III. If he judged there was any break in the profile, he pinpointed its radial and surface 
brightness coordinates in the /z-r profile. The characterization test was "blind", in the sense that the classifier did not 
know a priori the parameters that defined the models he was facing, and all models were "shuffled" , so as to avoid 
bias. The results of this simulation and the constraints they imply on the reliability of our results on analysis of real 
objects have been explained in section 3. 
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TABLE 1 

Results on the classification and characterization of Surface Brightness Profiles. 
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Note. — The name given for the objects is from GEMS (Rix et al. 2004). Coordinates are those retrieved by our software on 
GOODS-HST/ACS images (J2000.0). Rcdshifts arc from COMB017 (Wolf ct al. 2001, 2003). Stellar masses and B band absolute 
magnitudes arc extracted from Bardcn et al. (2005). The whole list can be consulted in the on-line edition. 
a The median value of the Pearson's R parameter for the linear fit from where hi is derived is fi — 0.993 
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TABLE 2 

Characterization of histograms in Figure 9. 



median 


(7 


z~0 23. 8f 


0.8 


O.K z < 0.5 22.3 


0.98 


0.5 < z < 0.8 21.5 


0.81 


0.8 < z < 1.1 20.6 


0.79 







Note. — Surface brightness at the break (psr ) m truncated galaxies, 
tg' band 



TABLE 3 

Characterization of histograms in Figure 10. 
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Note. — Rg r /hi ratio in truncated galaxies. 



TABLE 4 
Simulation Parameters. 
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Note. — Parameters for the creation of artificial galaxies (using GALFIT Peng et al. 2002), which enter our simulation on analysis of 
surface brightness profiles. The models were selected to have central surface brightness fi(r = 0) above 3cr a k y , and only models with Rb t I 
hi < 6, when appliable, are allowed. 

a hi >h 2 
h ht <h 2 
c R B r /hi < 6 



TABLE 5 

Results on Classification of Profiles. 
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a Numbcr of objects surveyed at each rcdshift range. 

^Number of objects analysed at each redshift range, after rejection of those not well suited to our study. 
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Fig. 1. — Example of the kind of break we define as "Jump Truncation", in which a minor discontinuity in intensity appears between the 
"inner" and "outer" exponential regions of the profile. In the left panel the discontinuous lines mark the position of the break according 
to the "Intersection Method" , while the continuous lines give the position given by the "Equal Deviation" method, specifically devised for 
these kind of profiles (16% of the objects which show a break). The ellipse on the image marks the position of the break, while the white 
line in the upper right corner is 1 arcsec long in the image. 
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Fig. 2. — Left: Absolute B-band magnitude of the galaxies in our sample against their estimated redshift. Right: stellar masses of the 
same objects against redshift. The vertical lines mark the limits of the 3 redshift bins in which we have divided our sample for study : 
0.1<z<0.5; 0.5<z<0.8; 0.8<z<l.l. The completeness levels in absolute magnitude and stellar mass at z~0.7 are -19.5 mag and 

M + ~ 3- 10 9 Mq respectively. 
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Fig. 3. — Comparison of the surface brightness profiles for 4 galaxies in common to the TP05 sample and ours. The vertical and horizontal 
lines mark the position of the Break, when there is such, according to each work (in the same color as in the profile). The differences 
between the two profiles as a function of the brightness meassured in TP05 are shown in the middle panels. We see how these differences 
grow sharply larger when we reach the ~26 mag/arcsec 2 isophote. On the right, images of each object are shown. Ellipses with the Break 
Radius as semi major axis are also drawn (horizontal white lines are 1" long). The ellipticities, position angles and centers of these ellipses 
are our own. 



18 



Azzollini, Trujillo & Beckman 



cn 

_o 

in 
Q 
O 

o 
o 



0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 



26 



. !. 



ct r = 0.2" 



-0.6-0.4-0.2 0.0 0.2 0.4 0.6 
UDF log(R Br TP (")) 



o 

CD 

if) 

u : 

^ 25 r 

cn 

D : 

Jl 24 r 



~ 3 - 23 r 
in 

Q : 

O ?? 
o 



r a, = 0.6 mag/arcsec z 



22 23 24 25 

UDF /j, B ™ (mag/arcsec 2 ) 



26 



Fig. 4. — Comparison of estimates of the Break Radius (Rbt \ left) and Surface Brightness at the Break (/iBr ; right) for a sample of 15 
galaxies (Type II), common to this work (ATB08) and to TP05. The discontinuous lines represent a 1:1 relation. 
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Fig. 5. — Results for simulated galaxies: left) classification success ratio and fraction of objects which are misclassifcd as being of the 
other 2 types as a function of the zg^g magnitude of the simulated galaxy; right) relative error in the visual estimate of iifl r , as a function 
of the input surface brightness at the "break" fj.Br in, for synthetic Types II and III (coded with different symbols). The uncertainty of 
the break position is ul' r = 5% for Type II objects, and (Tg^r = 13% for Type III. For both Types together, it is <j err = 8%. 
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Fig. 6. — Probability map of success in the classification of simulated galaxies, as a function of simulated hi and fiBr (up), and Rs r and 
[iBr (down). We present results for input models of Type II (left) and III (right). The points represent the position of real objects of our 
sample in the given planes. The 'crosses' mark the bins for which there were no input models available in the simulation. Note how the 
vast majority of Type II real galaxies are in a portion of the plane where the success classification is estimated to be very high (>80%). 
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Fig. 7. — Left: Break Radius [Rbt ) for "truncated" galaxies as a function of absolute B-band magnitude Mb , for 4 ranges of redshift. 
The data for z=0 are from PT06 (g'-band results). Right: Rb t at Mg =-20 mag, given by the best fit lines in the left panels, relative to 
z=0, for each redshift range. Ordinates are in logarithmic scale. A linear fit to log(i?5 r (z)/.Rg r (0)) - redshift gives a growth in Rb t by a 
factor 2.6±0.3 between z=l and z=0 for a same Mb ■ The discontinues straight line represents the relation for the local reference sample 
(the z=0 panel). The same applies to Figs. 8, 11, 12, 13, and 14. The numbers acompanying each point in the right panel give the 
population of objects which they represent. 
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Fig. 8. — Left: Break Radius (Rbt ) of "truncated" galaxies as a function of stellar mass M* , for 4 ranges of redshift. Local data are 
from PT06 (g'-band results). Right: R Br at M*=1O 1O M , given by the best fit lines in the left panel, relative to z=0, for each redshift 
range. Ordinates are in logarithmic scale. A linear fit to log(Rg r (z)/Rg r (0)) gives a growth in Rb t by a factor 1.3±0.1 between z=l and 
z=0, at a given M+ . The numbers acompanying each point in the right panel give the population of objects which they represent. 
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Fig. 9. — Evolution of the Surface Brightness at the Break (fj.Br )> in "truncated" galaxies, with redshift. Left: the fig r 's are rest- 
frame, i.e., they have been corrected for cosmological dimming (I(z) cx (l+z) — 4 ). From z~l to z~0 the median in the distributions of 
fj,g r has decreased by 3.3±0.2 mag/aresec 2 (a factor 20.9±4.2 in intensity). The median values are represented as vertical lines. For 
comparison we also show the median of the fiBr values given in Pohlen &; Trujillo (2006) for local galaxies, in g' band: fiBr 3 ,PT06 = 
23.8±0.8mag/arcsec 2 . Right: Median value of fiBr^or the distributions shown in the left panel, against their mean redshift. The larger 
error bars represent the standard deviation of the distributions, while the shorter ones give the error in the median values. The reported 
increase of 3.2±0.2 mag/arcsec 2 in /ig r with redshift is computed from the linear fit to the points. The numbers acompanying each point 
in the right panel give the population of objects which they represent. 



Cosmic Evolution of Stellar Disk Truncations 



21 



1 .2 



1 .0 

5 0.8 

cr 
CD 

£ 0.6 

CD 
> 

° 0.4 

CD 



0.2 



0.0 



0.1< = z<0.5 
0.5< = z<0.8 
0.8< = z<1 .1 
z~0 



2 3 

r b A 



Fig. 10. — Histogram of the ratio Rb t I hi for objects of the T-sample, in 3 different redshift bins, and the results in Pohlen & Trujillo 
(2006) at z~0 (g'-band results). The vertical lines mark the median values of Rbt / hi for each distribution. 
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Fig. 11. — Left: Rs r I hi for Type II galaxies as a function of absolute j3-band magnitude Mb , for 4 ranges of redshift. Local data from 
PT06 (cj'-band results). Right: Rsr I hi at Mb =-20 mag, given by the best fit lines in the left panels, relative to z=0, for each redshift 
range. The numbers acompanying each point in the right panel give the population of objects which they represent. 



Evolution at a given stellar mass 



_ z = 

. • "a 


. 0.1<z<0-5 


. 0.5<z<0,8 

,**••■ 

7 ' jf > . « j£ « - V— — — ~^ 


_ 0.8<z<1.1 









139 ~ - 


1 33 











0.6 



M. h 70 " 2 [10'°M B ] 



M. h 70 "'[10'°M B ] 



Fig. 12. — Left: Rs r I hi of "truncated" galaxies as a function of stellar mass M* , for 4 ranges of redshift. Local data are from PT06 
(g'-band results). Right: Rbt I hi at M* =10 10 Mq , given by the best fit lines in the left panel, relative to z=0, for each redshift range. 
The numbers acompanying each point in the right panel give the population of objects which they represent. 




Fig. 13. — Left: hi for Type II galaxies as a function of absolute B-band magnitude Mg , for different ranges of redshift. Local data from 
PT06 ((/-band results). Right: hi at Mg =-20 mag, given by the best fit lines in the left panels, relative to z=0, for each redshift range. 
Ordinates are in logarithmic scale. A linear fit to log(iJg r (z)/Rg r (0)) - redshift gives a growth by a factor 1.4±0.5 in hi between z=l and 
z=0 at Mg =-20 mag, a result compatible with no evolution. The slash-dotted line is the relation between size h (equivalent exponential 
scale length deduced from the Sersic effective radius in r' band) and absolute magnitude found by Shen et al. (2003) in local late-type 
(Sersic n<2.5) objects using the SDSS. Note how, as expected, hi is larger than h for all the galaxies. The numbers acompanying each 
point in the right panel give the population of objects which they represent. 
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Fig. 14. — Left: hi of Type II galaxies as a function of stellar mass M* , for 4 ranges of redshift. The data for z=0 is from PT06 (g'-band 
results). Right: hi at M+ =10 Mq , given by the best fit lines in panel on the left, relative to z=0, for each redshift range. Ordinates are in 
logarithmic scale. A linear fit to log(Rg r (z)/Rg r (0)) gives a shrink by a factor 0.8±0.1 in hi between z=l and z=0 for a same M* . The 
slash-dotted line is the relation between size (equivalent exponential scale length) and stellar mass found by Shen et al. (2003) for local, 
late-type (Sersic n<2.5) objects in the SDSS (taken from Fig. 11 in their paper). Note how hi is larger than h for all the galaxies. The 
numbers acompanying each point in the right panel give the population of objects which they represent. 
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